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Abstract 

The  dispersion  spectrum  is  found  for  axially  symmetric  torsional  waves 
in  an  elastic,  bimetallic  rod  with  cylindrical  core  and  concentric  outer 
casing.  Plotting  of  the  various  branches  of  the  spectrum  is  simplified  by 
the  presence  of  discrete  invariant  points  which  are  independent  of  the 
material  properties,  and  through  which  the  spectral  lines  must  pass.  The 
slope  and  curvature  of  the  spectral  lines  at  cut-off  frequencies,  the 
asymptotic  approximations  at  high  frequency,  the  non-existence  of  complex 
branches,  the  problem  of  co-existence  and  the  concept  of  energy  flow  are 
also  studied. 
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Introduction 

I 

The  problem  of  wave  propagation  in  a  piecewise,  homogeneous  elastic 

i 

solid  is  of  interest  in  a  number  of  physical  areas,  e.g.,  [1,2,3].  This 
study  concentrates  on  a  geometrically  simple  problem,  which  contains  many  i 
aspects  of  more  general  problems  but  has  not  been  fully  explored  in  the 
literature.  In  particular,  the  dispersion  equation  is  obtained  and  analyzed 
for  axially  symmetric  torsional  waves  in  a  bimetallic  circular  cylinder 
with  solid  cylindrical  core  perfectly  bonded  to  a  cylindrical  casing  of  a 
different  elastic  material. 

Earlier  investigations  in  this  area  are  due  to  Armenakas  [4],  Reuter 
[5],  and  Haines  and  Lee  [6].  Using  Pochhammer's  equation  for  torsional 
waveB  in  elastic  cylinders  [7],  Armenakas  presented  numerical  results  for 
a  few  typical  choices  of  material  and  geometric  parameters  for  real  wave 
numbers  only.  Reuter  carried  out  an  asymptotic  analysis  for  the  low 
frequency  long  wave  length  region  and  determined  the  expression  for  the 
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phase  velocity  of  the  lowest  torsional  mode.  Haines  and  Lee  reinvestiga¬ 
ted  the  problem  and  numerically  demonstrated  the  properties  of  the  dis¬ 
persion  spectrum  for  a  general  composite  cylinder  for  both  real  and 
imaginary  wave  numbers.  They  also  formally  pointed  out  the  absence  of 
complex  segments  of  the  dispersion  spectrum  and  rederived  the  low  fre¬ 
quency  long  wave  length  approximation  for  the  lowest  torsional  mode. 

The  present  study  amplifies  and  further  extends  the  research  carried 
out  by  these  earlier  investigators  in  this  problem  area.  In  addition  to 
a  solution  for  the  actual  frequency  equation,  obtained  by  numerical  solu¬ 
tion  of  transcendental  eigenvalue  equation,  a  number  of  other  significant 
results  have  been  obtained.  First,  it  is  shown  that  the  eigenvalue 
equation  is  separable  in  the  Bernoulli  sense  so  that  it  may  be  written  as 
the  sum  of  two  analytic  functions,  one  of  which  is  a  function  of  the  radial 
casing  wave  number  and  the  other  is  a  function  of  the  radial  core  wave 
number.  The  separation  constant  Q  ,  represents  the  elastic  coupling  at 
the  bonded  Interface  between  the  core  and  casing  and  therefore  can  be 
incorporated  in  the  Interface  continuity  equations.  The  solutions  to  the 
eigenvalue  problem  is  thus  an  analytic  function  of  the  Lie-parameter  Q  . 

In  particular,  the  lower  limiting  value  Q  -*■  0+  corresponds  to  a  traction- 
free  boundary  eigenvalue  problem  and  the  upper  limiting  value  Q  +  • 
corresponds  to  a  displacement- free  boundary  eigenvalue  problem.  The 
problem  of  perfect  bonding  at  the  interface  is  therefore  an  intermediate 
eigenvalue  problem.  It  therefore  follows  from  Weyl-Courant  minimax  theorem 
in  the  theory  of  self-adjoint  differential  operators  that  the  Rayleigh 
quotient  is  a  continuously  increasing  function  of  the  Lie-parameter  Q  , 
and  therefore  the  eigenvalues  of  the  intermediate  problem  are  nested  by 
the  eigenvalues  of  the  two  limiting  cases,  that  is,  (*)_«+  <  u  < 
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The  dispersion  equations  for  these  limiting  rases  take  on  a  verv  simnlr 
form  and  can  be  easily  plotted  and  there! ore  serve  as  bounds  on  the  dis¬ 
persion  spectrum  of  the  intermediate  problem,  as  required  by  the  minimax 
theorem. 

To  qualitatively  sketch  the  dispersion  spectrum,  we  have  also  deter¬ 
mined  analytically  the  cut-off  frequencies,  their  slope,  and  their  curvature 
at  long  wave  length.  We  have  also  shown  that  except  for  the  lowest 
torsional  mode,  coincidence  of  frequencies  at  cut-off  does  not  exist  and 
therefore,  in  general,  all  branches  except  the  lowest,  have  zero  slope  at 
infinite  wave  length.  The  lowest  torsional  mode  is,  however,  an  exception 
because  at  cut-off,  zero  frequency  is  a  non-simple  double  root  of  the 
eigenvalue  equation  and  two  linearly  independent  eigenfunctions  can 
co-exist  at  this  frequency,  giving  rise  to  coincidence  phenomena,  and  non¬ 
zero  slope.  We  have  further  shown  that  in  the  long  wave  length  region,  the 
curvature  of  the  branch  is  greater  (less)  than  c|/w  when 
greater  (less)  than  '  ^us  ^or  a  typical  set  of  material  parameters, 

all  branches  of  the  spectrum  except  the  lowest,  have  similar  geometric  pro¬ 
perties  in  the  neighborhood  of  cut-off  frequencies.  To  obviate  elaborate 
computations  in  the  high  frequency  region  we  have  also  obtained  McMahon- 
asymptotic  representations  of  the  various  frequency  equations  which  enter 
in  the  discussion.  We  finally  conclude  this  study  with  a  brief  discussion 
of  group  velocity  and  flow  of  energy  flux  in  bimetallic  cylinder. 
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Basic  Equations 

For  torsional  waves  in  a  homogeneous,  isotropic,  elastic  medium 
the  displacement  equation  of  motion  in  an  open  domain  R  is  given  by 


1 

77 


u  + 


a2u  i  a2u 

177  -77  377  • 


c2  >  0 


(1) 


where  c 2  *  y/p  is  the  square  of  wave-speed  of  the  shear  waves  in  an 
Infinite,  homogeneous ,  isotropic  elastic  medium  with  rigidity  modulus 
p  and  mass  density  p  .  The  tangential  displacement  component  is 
u  =  u0(r,z;t),  where  r  is  the  radial  direction,  z  the  longitudinal 
direction  and  t  is  the  time  parameter,  e.g.,  [7].  Equation  (1)  is 
a  two-dimensional  wave  equation  and  admits  plane  waves  propagating  in 
the  longitudinal  z-direction.  Assuming  solutions  of  the  form 


u(r,z;t)  -  V(r)  exp  i(£z-wt) ,  (i  = 


(2) 


we  find  that  V(r)  must  satisfy  the  equation 


1  d_ 
r  dr 


(r  -  iy  V  +  <2V  -  0  , 


K2  $  0 


(3) 


where 


k2  =  (w/c)2  -  C2,  (4) 

* 

is  the  square  of  radial  Wave  number,  £  is  the  wave  number  in  the 
longitudinal  direction,  u>  is  the  angular  frequency  with  real  period 
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2n/w  and  frequency  u/2ir  .  The  radial  function  V(r)  is  the  solution 
of  Bessel  equation  of  order  one  when  k  i  0  ,  and  in  the  limiting 
case  when  k  ■  0  ,  the  Bessel  equation  reduces  to  an  Euler  equation. 

In  order  to  accommodate  the  limiting  case  when  the  radial  wave  number 
approaches  the  value  zero,  we  take  the  general  solution  of  Eq.  (1)  in 
the  form 

ue  ■  [A  7  -  f 

where  A  and  B  are  suitable  coefficients  to  be  determined  from  the 
boundary  conditions. 

We  define  our  bimetallic  composite  cylinder  R:F:,^  x  R2  to  be  the 
union  of  two  cylinders  R1  and  R2  such  that 

Rj  :  0  <  r  <  ro  ,  |*|<- 

R2:ro<r-a*  M<“ 


] 


i  exp  i(£z-wt) 


(5) 


a  >  rQ  with  interface  r  ■  r  ,  which  is  assumed  to  be  perfectly  bonded. 
The  material  properties  of  the  two  cylinders  are  assumed  to  be 


R1  :  ^1*^1^  »  R2  :  ^2*p2^  * 


The  shear  stresses  are  given  by 


Tze  m  u  3z  ’ 


3  rl  \ 

\emvr  3 7(7  ue)  * 


(6) 
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and  therefore  using  Eq.  (5)  we  find  that 


Tze  ■  ivS|A~  ^Ur)  -  f  icBY^(tcr)jexp  i(£z-u)t)  , 
Tr0  "  4a  J2^Kr^  ~  f  K2BY2(Kr)lexp  i(Cz-ut)  . 


(7) 


For  a  solid  cylinder  in  region  R^  ,  the  boundary  conditions  as 
r  -*•  0+  require  that  lu^1^  (0+,z;t)  |  and  I t (0+,z;t)  |  be  bounded 
and  therefore  suitable  forms  for  displacement  and  shear  stress  are 


u<»  -  —  A1J1(<1r)  exp  i(Cz-(dt)  , 
Tr0)  "  ~2^lAlJ2^Klr^  exp  , 

where 


(8) 


<1  =  [(m/Cj)2  -  S2]*2  , 


(9) 


C2  5  Vj/Pj.  • 


For  a  hollow  cylinder  in  region  R2  ,  the  suitable  forms  for  dis¬ 
placement  and  shear  stress  are 


.(2) 


(*2r)  -  j  k2B2Y 


l(K2r)] 


exp  i(£z-«t)  , 


-  -y2|"2A2J2(tc2r)  -  K2B2Y2(ic2r)j  exp  i(Cz-ut)  , 


(10) 


where 


<2  =  [(w/c2)2  -  C2]  ^2, 
c2  =  V2/p2  * 


(ID 


Frequency  Equation 

For  a  bimetallic  cylinder  with  generators  parallel  to  the  z-axis 
the  boundary,  interface  and  boundedness  conditions  are 


1.  T 


2. 


u) 

1  m 

0 

3R2 

re 

9R2 

(1) 

x  (2) 

re 

3R1 

re 

3R2 

3R1 

at 

(2) 

ue 

3R2 

3R1 

3R2 

3R,  :  (r  ■  a)  x  z  , 


:  (r  -  rQ)  x  z  , 


(12) 


3. 


u<*> 

,  <  M  , 

x(l) 

e 

+  * 

r-K) 

re 

<  M 


r-K) 


We  assume  that  R^  and  R^  are  perfectly  bonded  at  their  common 
interface  R^Or2  :  (r  “  tQ)x  z*  therefore  assume  that  the  ensuing 
motion  of  the  bimetallic  cylinder  will  have  a  common  frequency  w/ 2ir 
and  a  common  spatial  period  2ir/£  in  the  axial  direction.  However, 
the  radial  wave  numbers  e  R^  and  k2  e  R2  will  in  general  be 
different.  The  two  wave  numbers  and  <2  are  not  independent  and 


only  those  values  are  admissible  for  which  the  boundary  and  Interface 


conditions  (12)  are  satisfied.  We  have  one  boundary  condition  (12)^ 
and  two  continuity  conditions  (12)2>  From  Eqs.  (8) ,(10)  and  (12)  we 
find  that  there  are  five  unknowns  and  <2  •  Since  we 

have  only  three  homogeneous  equations  relating  the  five  unknowns,  it 
follows  that  the  two  radial  wave  numbers  and  <2  are  not  indepen¬ 

dent.  For  the  three  homogeneous  equations  to  have  a  non-trivial  solution 
the  determinant  of  the  coefficients  A^,A2  and  B2  must  be  zero.  This 
leads  to  the  frequency  equation 

na2P2  -  oy[J2(Y)/J1(Y)3(aQo+  2Sq)«  0  ,  (13) 


where 


Y  =  r  k 


o  1 


T1  =  rQ/a  , 


a  =  2  » 


and 


P-  =  P  +  —  Q  +  —  R  +  —V  S 
2  o  no  o  a  o  na'1  o 


The  functions  P  ,P~;  Q  ,R  and  S  are  defined  as 
0200  o 


Pn  5  “  Va)Jn(ria) 

<?0  5  J0(a>Yo(na)  -  Jg(no)Y0(a) 
Ro  =  j;(a)YQ(na)  -  J0(no)Y;(o) 
SQ  =  j;(a)Y;(na)  -  j;(na)Y;(a) 


n  -  0,2 


» 


(14) 


» 


(15) 
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where  primes  indicates  differentiation  with  respect  to  their  respec¬ 
tive  arguments,  e.g.  [8], 

When  a  is  imaginary  and  y  real,  i.e.,  when  a  -*•  ia  and 
y  -*■  y  ,  the  frequency  equation  (13)  takes  the  form 

T)a2?2(a)  -  oylJ2(y)/J1(Y)](aQo(a)  +  2SQ(a))  =  0  ,  (16) 

where  now  a  and  y  are  both  real  and 

P2(o)  =  I2(a)K2(no)  -  I2(na)K2(o)  , 

Qo(a)  =  IQ(a)K1(na)  +  I1(no)K0(o)  , 

s  r (n«)Ko(a)  -  I0(o)K’(na)  ,  (17) 

S0(a)  2  I1(na)K1(a)  -  I1(a)K1(r»ot)  , 

=  r(a)K;(na)  -  i;(no)K;(a)  . 

When  both  a  and  y  are  imaginary,  i.e.,  a  ■+■  ia  ,  y  iy  the 
frequency  equation  (13)  takes  the  form 

na2P2(a)  +  oy[I2(y) /ij^Cy)  ]  (aQQ(o)  +  2SQ(a))  =  0  ,  (18) 

where  a  and  y  are  both  real  in  Eq.  (18) ,  I  and  K  are  modified 

n  n 

Bessel  functions  and  other  symbols  are  same  as  defined  earlier. 

The  case  of  grazing  incidence  occurs  when  either  y  =  0  or  when 
a  ■  0  .  When  a  -*■  0  ,  the  frequency  equation  takes  the  form 

yJ2(y)  -  0  . 


(19) 


When  y  Q  ,  the  frequency  equation  takes  the  form 


aP2(a)  "  0  • 

Equation  (20)  has  no  real  zeros  except  a  =  0,  since 

P2(a)  =  I2(a)K2(na)  -  I2(na)K2(a)  , 

•*»  2  cosh  a(l-n)+  •••  ,  a  >>  1  , 


(20) 


(21) 


which  has  no  zero>.  Hence  the  dispersion  spectrum  never  intersects 
the  line  y  ■  0  ,  but  approaches  it  asymptotically  from  above.  This 
eliminates  the  existence  of  imaginary  values  of  a  and  y  in  Eq.  (13). 
However,  Eq.  (19)  admits  simple  zeros  and  therefore  in  the  region 
w/c2  <  ?  <  w/c^  ,  y  is  real  and  a  is  imaginary. 

The  last  case  is  when  a  is  real  and  y  is  imaginary,  i.e.,  when 
a  -*■  a  and  y  ir  .  The  frequency  equation  in  this  case  takes  the  form 

na2P2(a)  -  aY[I2(Y)/l1(Y)](aQ0(a)  +  2SQ(a))  »  0  ,  (22) 

where  now  a  and  y  are  both  real,  and  other  symbols  are  as  defined 
earlier.  This  equation  replaces  Eq.  (16)  when  the  wave  speed  in  the 
core  is  greater  than  that  in  the  casing.  Such  a  situation  is  certainly 
possible,  but  is  of  less  physical  interest. 


Analysis  of  Frequency  Equation 

The  frequency  equation  (13)  can  be  treated  as  a  function  of  two 
independent  variables  a  and  y  and  has  an  interesting  property  that 
the  equation  is  separable.  It  is  easy  to  see  that  this  equation  can  be 
written  in  the  form 

3t(o)  -  yJ2(y)/j1(y)  ,  (23) 

where 

na2P2(a) 

y (a)  H  a[aQo(a)+2So(oO]  * 

Since 


£*»>  -°  • 

(24)* 

^•(yJ2(y)/J1(y))  -  o  , 

it  follows  that  there  exists  a  separation  constant  Q  such  that 


^(a)-  Q  a  0  » 

YJ2(y)  -  QJL(Y)  -  0  , 


(25) 


where  0  <.  Q  <  <*>  and  can  be  considered  as  a  continuous  parameter  e.g. 

[9,10]. 


*  Note  that  this  process  of  separation  is  capable  of  generalization. 
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The  second  of  Eq.  (25)  corresponds  to  the  case  of  a  solid  cylin¬ 
der  with  elastic  boundary  condition  at  r  =  rQ 

■  -OV'o  •  (26) 

When  Q  -*■  0+,  we  have  the  case  of  a  traction-free  solid  cylinder  and 
q  +  w  corresponds  to  the  case  of  a  solid  cylinder  with  rigid  encastered 
surface.  The  rigid  boundary  conditions  slowly  change  to  traction-free 
boundary  conditions  as  we  gradually  change  the  parameter  Q  from  <» 
to  0+  . 

Similarly  the  first  of  Eq.  (25)  corresponds  to  the  case  of  an 
hollow  cylinder  with  elastic  boundary  conditions 


i>  T<2)/u<2)  -  -0Ul/ro,  (r  -  r0>  .  . 

ii)  Tr0)  “  0  (r  -  a)  x  z  . 


(27) 


By  varying  Q  from  0+  to  <*>,  the  interface  boundary  condition  at 
r  ■  rQ  changes  from  traction-free  to  rigid  condition.  Thus  when  Q  ■  0, 
we  get  from  Eq.  (25) ^ 


P2(a)  -  0  if  a  i4  0  , 


(28) 


which  is  the  frequency  equation  of  a  traction-free  hollow  cylinder. 
When  (J  f  *  ,  Eq.  (25) ^  gives  us  the  frequency  equation 


aqo(o)+2S0(a)  -  0 


if  o  0  , 


(29) 


which  is  the  frequency  equation  of  a  hollow  cylinder  with  traction- 
free  external  surface  and  rigid  internal  surface. 


We  therefore  conclude  that  when 
tions  are 


0+, 


the  two  frequency  equa- 


i)  J2(y)  -  0  Y  ¥  0  , 

ii)  P2(a)  -  0  a  4  0  , 

and  when  Q  f  «°  ,  the  two  frequency  equations  are 


(30) 


i)  Jj/y)  -  0  , 

ii)  aQo(a)  +  2S0(a)  *  0  if  o  +  0. 


(31) 


It  can  be  easily  verified  that  these  frequency  equations  have  only  simple 
zeros.  Furthermore,  these  are  also  the  zeros  of  the  general  frequency 
equation  (13).  In  addition  these  simple  zeros  have  the  interesting 
property  that  they  do  not  depend  upon  the  choice  of  the  rigidity  ratio  o  . 
Thus  all  of  these  roots  of  the  frequency  equation  are  invariant  with 
respect  to  the  ratio  of  the  rigidity  moduli!  of  the  bimetallic  cylinder. 

Let  Yp  i  P  *  1,2,3,...  and  ,  q  «  1,2,3,...  be  the  roots 
of  the  two  frequency  equations  (30).  Then  using  these  values  of  a  and 
Y  ,  the  spectral  branches  of  the  two  frequency  equations  are  given  by 
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(aw/c2)^  -  a2  +  (a?)2  , 
(au>/c2)2q  -  |-[(Yp/n)2  +  (a£)2]  , 


(32) 


where  p  = 


(aC) 


2 

pq 


p^/p^  .  The  points  of  intersection  of  these  branches  are 

o(rp/n)2  -  pa2 
p  -  a  * 

o[(Yp/n)2  -  a2] 
p  -  o 


(33) 


From  the  earlier  discussion  it  is  obvious  that  these  points  of  inter¬ 
section  are  invariant  points  of  the  dispersion  spectrum  through  which 
the  branches  of  the  frequency  equation  (13)  must  pass.  Furthermore, 
these  are  the  only  points  at  which  these  branches  can  pass  through 
these  Q  -»■  0+  and  Q  +  “  spectral  branches  which  may  then  be  thought 
of  as  bounding  curves. 

For  large  values  of  the  argument,  frequency  equations  (30)  have  the 
Hankel-type  asymptotic  form 


i)  tan(y-ir/4)  ■  y  , 

ii)  tana(l-n)  "  p-  . 

o  qa 


(34) 


It  is  relatively  easy  to  find  the  roots  of  these  transcendental  equa¬ 
tions.  These  roots  when  used  in  Eq.  (33)  determine  the  invariant  inter¬ 
section  points  of  the  frequency  equation  (13)  in  the  high  frequency 
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Now  consider  the  remaining  two  frequency  equations  (31) .  For 
large  values  of  the  argument,  the  Hankel-asymptotlc  forms  for  these 
equations  are 


i)  tan(Y-Tr/4)  +  —  -  0  , 

ii)  tana(l-n)  +  "  0  • 


(35) 


Let  Yn  ,  n  ■  1,2,3,...  and  am  *  m  *  1,2,3,...  be  the  zeros  of 

these  two  frequency  equations  which  for  large  values  of  y  and  a 

n  m 

agree  with  the  zeros  of  Eq.  (31) .  Then  in  the  high  frequency  region 
the  spectral  branches  of  the  two  frequency  equations  (31)  are  given 
by  Eq.  (32)  if  we  replace  y^  by  Yn  and  by  aR  .  Under  the 

same  transformation,  the  points  of  Intersection  of  the  branches  are 
given  by  Eq.  (33).  Again,  these  additional  intersection  points  are 
invariant  points  of  the  spectrum  through  which  the  branches  of  the 
frequency  equation  (13)  must  pass. 

Now  Q  is  a  continuous  parameter  and  for  various  values  of  Q 
one  may  easily  find  the  roots  of  the  frequency  equation  (25) .  Because 
the  frequency  equations  are  entire  functions,  it  is  possible  to 
express  a  and  y  in  terms  of  power  series  in  Q  and  one  can  then 
study  the  properties  of  the  a  and  y  trajectories  with  respect 
to  Q  .  However,  such  an  analysis  is  not  of  much  interest  because 
the  roots  a  of  Eq.  (25)  ^  depend  upon  the  choice  of  ratio  o  and  the 


intersection  points  of  the  a-  and  y-  branches  are  not  invariant 
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points  of  ths  spectrum.  With  s  viov  towards  ussfulnsss  in  ths  study  of 
Floqust  waves  in  blaetallic  cylinder  with  periodic  structure,  only 
these  invariant  points  narlt  special  consideration.  As  such, no  attempt 
is  made  to  study  the  properties  in  the  neighborhood  of  the  o-varlant 
intersection  points  for  other  values  of  parameter  Q  . 


1 
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Long  Wavelength  Properties 

With  the  invariant  points  at  our  disposal,  the  branches  of  the 
spectrum  pertaining  to  the  frequency  equation  (13)  can  be  qualitatively 
traced  if  one  has 

a)  the  cut-off  frequencies, 

b)  the  slope  of  the  branches  at  cut-off,  and 

c)  the  curvature  of  the  branches  at  cut-off. 

Supplementing  this  information  with  the  properties  at  long  wavelength 
and  low  frequency;  zeros  in  the  zero-frequency  plane;  and  the  zeros  at 
grazing  incidence,  provides  sufficient  information  to  trace  the  branches 
of  the  spectrum  qualitatively.  Additional  information  based  on  the  theory 
of  differential  equations  reveals  that  the  problem  of  coexistence  does 
not  exist  because  all  the  roots  are  simple  and  furthermore,  complex 

roots  are  inadmissible.  However,  the  lowest  torsional  mode  is  an  exception. 


i)  Slope  at  cut-off 

If  we  represent  symbolically  the  frequency  equation  (13)  by  the 
functional  form  F(a,y)  ■  0,  then  using  the  differentiation  formulas  of 
implicit  functions,  we  get 


&oi 


F  a.  +  F  Yr 

_o_5 _ X_L 

Fa  +  F  y 
a  u  Y  u) 


(36) 


where  a  =  a(u>(C)»0,  Y  =  y(ui(C)»C)  and  F^.F^.o^ , . . .  indicate  par¬ 
tial  derivatives  with  respect  to  the  variables  and  u>  , 

respectively.  How 
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<*5  •  -  *2C/a  , 


(37) 


and  In  the  limit  as  £  0,  a^  |  -*•  0  and  |  q  -*■  0  .  Therefore  it 
follows  that  at  cut-off  frequencies,  the  slope  of  the  branches  Is  in 
general  zero,  i.e.. 


5*0 


0  , 


(38) 


except  when  the  denominator  (Fa  +  F  y  )  -  0.  For  non-zero  values  of 

aw  y  w 

wave  number  £  the  slope  of  the  branches  will  vanish  if 


1  3F  .  2  1  3F 

a  3a  n  y  3y 


0  . 


(39) 


In  the  region  of  long  wavelength  and  low  frequency,  when  both  £  and 
w  tend  to  zero,  Eq.  (36)  takes  an  indeterminate  form,  since  it  can 
be  rewritten  as 


3w 

at 


YFa  +  n2aFY  ,  £ 

■yc1Fo  +  n‘c2aV  "  ’ 


(40) 


This  can  be  written  as 

s»2  .  T?«  1 

a?3’  yc1f0  +  n*c2aI,7  ’ 

which  suggests  that  in  the  low  frequency  long  wave  length  region,  w2  is 
linearly  proportional  to  £2.  To  determine  the  slope  of  the  branch  in 
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this  region,  when  both  w  <<  1  and  £  <<  1  ,  we  use  power  series 
expansion  of  Bessel  functions  for  fixed  order  and  small  values  of  the 
argument .  The  frequency  equation  In  this  region  takes  the  approximate 
form 


any2  +  (rT1-  n3)a2  *  0  , 


or  equivalently 


_1£i£2ILL 

l-(l-p)n4 


(41) 


as  conjectured  earlier.  From  here  we  immediately  conclude  that  the  slope 
of  the  lowest  branch  at  long  wave  length  and  low  frequency  is  given  by 


du 

fi  ,  (°-p)n4 

3£ 

>  "C2 

[  +  l-(l-p)n4J 

o 


_  ,  w  «  1 

V  £  «  1 


(42) 


which  agrees  with  the  limiting  value  obtained  by  Reuter  [  5  ]. 

From  Eq.  (36)  we  also  see  that  the  slope  is  indeterminate  when 


F  a.  +  F  y_  ■  F_ 
a?  Y  5  £ 

Fa  +  F  y  ■  F 
a  u  y  w  w 


0  . 
0  . 


Simultaneous  solution  of  the  three  equations 


F  -  0  , 


F_  »  0  and  F  -  0  , 
£  u 


(43) 
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determine  u  ,  5  and  o  for  which  we  have  coalescence  of  spectral  lines. 
These  are  points  of.  coincidence  and  as  shown  in  Ref.  [Ill  the  slope  at 
these  points  is  given  by 


9ai  , 
«5  35 


(44) 


which  gives  us  the  value  of  slopes 

(•“)F  *  (-*  r  ±^2r  ~  F  Frr>M  •  (45) 

l3v  wo)  [  0)C  oj£  oko  J 

When  F  I  -*■  0  ,  as  in  the  present  case,  the  slope  of  the  branches  at 
points  of  coincidence  takes  the  simple  form 


In  the  present  problem  it  can  be  verified  that  Eq.  (43)  has  only  one 
non-trivial  solution  (01,4) : (0,0)  ,  which  reveals  that  the  origin  in  the 
(a), 4) -plane  is  a  zero  of  multiplicity  2  .  In  the  6-neighborhood  of  the 
origin  F^  2  ,  ■  -c|[l  +  (o-p)n‘*/[l  -  (1— p)  13** ]  1 »  and  therefore  from 

Eq.  (46)  we  get  (3u/3 O^0  Q)»  which  agrees  with  the  result  derived  earlier 
for  the  lowest  branch.  As  remarked  earlier,  this  multiplicity  leads  to 
co-existence  of  two  linearly  independent  eigenfunctions. 


ii)  Curvature  at  cut-off 

At  cut-off  frequencies  when  the  slope  is  zero,  the  curvature  of 
the  branches  is  given  by 


w 


F  att  +  F  y„_] 
a  M  Yt€5 

4-*o 

Fa  +  F  v 
a  w  ytw 

Kr*> 

«D  -M> 


(47) 


J 

5 
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a  -  -  -ac2/w  , 

^  5-*o 


r«  „  ’  •r»ci/"' 


x  -  a/c0  » 

0)  „  Z 

w  C->o 


*.  '  r=/Cl  ’ 


and  therefore  the  curvature  of  all  those  branches  at  cut-off,  whose 
slope  is  zero,  is  given  by 


C2  ^(C-./Ca"C*/C-) 

“  1  +  Fa/F  +n(c2/c'“)  * 


For  high  frequencies  and  long  wavelengths,  both  a  and  y  are 
large.  Using  asymptotic  expansion  of  Bessel  functions  of  large  argument 
and  fixed  order,  the  frequency  equation  (13)  takes  the  Hankel-asymptotic 


nasina(l-n) (cosy-siny)  +  aycoso(l-n) (cosy  +  siny)  s  0  . 


From  here  one  can  easily  show  that 


yj5-K) 


;2y2+  n2q2)+(a2Y2-  n2a2)sin2 
naa(2y  +  cos2y) 


f-noycosi 


and 


(1-n)  ~{ (l+op)+(l-op)sin2Y  }+  2y  Sap 
n  2  o  o 


p  c  V riy  —  l  V -L-ru 

-Z  +  n  -1  *  a  _ 2 _ 

Jy  clJ^  o  (2yQ  +  cos2y  ) 


Therefore  the  curvature  of  the  higher  branches  at  cut-off  is  given  by 


g2u  c2  •  (o/p)  ^2(o-p)(2yo+  cos2yo) 

C**°  W  2y  (op)  ^  +(l-n)  (~)  { (l+op)+(l-op) eln2y 


,  (50) 


where  yo  =  n(p/o)  ^(au/c2),  and  ao>/c2  is  a  dimensionless  frequency. 


(l+op)+(l-op)sin2Yo  ■  (cosY0+sinyo) 2+  op (cosYo-sinYo) 2  , 


and  is  therefore  a  positive-definite  quantity.  Also,  for  y  >~  and 

o  2 

real,  2yo  +  cos2yq  =  (2yq+1) cos2Y0+  (2yq-1) sin2YQ  >  0  .  Hence  the 
second  term  in  the  parentheses  is  positive  if  a  >  p  and  negative  if 
o  <  p  .  When  o  ■  p  ,  the  curvature  is  c^/w  ,  which  is  the  curvature 
of  the  a-branch  at  cut-off  frequency.  Thus  starting  at  cut-off  fre¬ 
quency  with  zero  slope,  the  branch  of  the  spectrum  lies  above  (below)  the 
a-hyperbola,  when  a  is  greater  (less)  than  p  . 


iii)  Cut-off  frequencies 

To  sketch  the  higher  branches  qualitatively,  it  now  remains  to 
determine  the  cut-off  frequencies.  At  £  “  0  ,  these  frequencies  are 
given  by  the  higher  roots  of  the  frequency  equation 


tan(l-ri)n  -(op) ^2tanjji(p/o)  “  f]  " 


where  Q  =  am/c2  -  c|^0  -  ± 


Asymptotic  Formula  for  the  Zeros  of  Frequency  Equation  when  j Q |  -*■  “ 

The  zeros  of  the  frequency  equation  (30)  when  Q  =  0  are  well  tab¬ 
ulated  and  their  asymptotic  form  can  also  be  found  in  [8].  Also,  one 
can  find  the  zeros  of  J^(y)  *  0  ,  well  tabulated.  However  asymptotic 
formulas  for  the  zeros  of  frequency  equation  (31) ^  are  not  available  and 
in  this  section  we  obtain  their  McMahon-asymptotic  representation  p .2], 
Consider  the  equation 

aQo(a)  +  2SQ(a)  -  0  ,  (52) 


where 


Q0(«)  “  J1(na)YQ(a)  -  JQ(a)Y1(na)  , 

SQ(a)  »  Y1(na)J1(a)  -  Y^a)  J^na)  , 

and  n  =  rQ/a  14  ^  •  Using  Hankel -asymptotic  expansion  of  cylinder  func¬ 
tions  for  fixed  order  and  large  arguments,  it  can  be  shown  that  Eq.  (52) 
can  be  written  as 

cos[(l-r))a  +0]  -  0  ,  (53) 


where 


c 080  'u  1  +  (1-  +  6f|  -  7rj2)  +  ...  , 

8lne  '  sib  +  TTsibr  °  +  15n  +  2V  +  21"3)  +  •••]• 


(54) 


The  zeros  of  Eq.  (53)  are  given  by 


(l-n)o  “  8  -  0,  (55) 

n  n 

where 

8  -  (2n-l)  |  •  n  =  1,2,3,...  (56) 

n  z 

and  6  is  the  solution  of 

tan6  *  8fe[(5n"1)  +  78^)2  <5  +  9ri  "  A5ti2  +  63r,3)  +  ’  ‘  *]  *  (57) 

When  a  »  1,  it  follows  that  tan6  <<  1  and  therefore  we  can  write 

0  ■  tan0  -  — ■  tan30  +  tan50  -  ...  (58) 

From  Eqs.  (57)  and  (58)  we  therefore  get 

e  ■  sfe  +  T&3  (7  -  15,l3)  +  •  •  •  •  <59> 

From  Eqa.  (55)  and  (59)  it  now  follows  that  for  large  values  of  a,  the 
asymptotic  value  of  the  roots  is  governed  by  the  algebraic  equation 

a  ■  y  +  p/a  +  q/a3  +  ...  ,  (60) 


where 


Using  Lagrange's  expansion  theorem,  the  roots  of  the  algebraic  equation 
(60)  are 

a  Y  +  p/Y  +  (q-P2)/Y3+  •  (62) 


Hence  for  large  values  of  a  ,  the  asymptotic  formula  for  the  zeros  of 
the  frequency  equation  (52)  have  the  explicit  representation 


3(l-5n)  _  12(7-n-60n2+  135n3+  I5n4) 

“n  Yn  8riyn(l-n)  (8nyn) 3  (1-h) 2 


0(1/Y3) 


(63) 


where 


26 


Inadmissibility  of  Complex  Wave  number  £ 

In  this  section  we  show  that  in  the  case  of  a  bimetallic  rod,  the 
wa^e  number  £  can  be  real  or  imaginary,  but  cannot  be  complex.  This 
easily  follows  by  using  the  theory  of  singular  Sturm-Liouville  eigen¬ 
value  problems  [13]  . 

Consider  Eq.  (3)  which  can  be  written  in  the  self-adjoint  form 
Lu  +  «2ru  *  0  ,  r  e  R  :  R^x  Rj,  k2  ^  0  (65) 

where 


and  r  “  0  is  a  singular  point.  The  domain  of  the  operator  L  is 
the  class  of  C2(I)  functions  which  satisfy  the  boundary  and  continuity 
conditions  (12) . 

Let  u^(r)  and  ^(r)  be  two  complex-valued  functions  which 

satisfy  the  equations 

* 

i 

| 

'  Lu.  +  <?ru.  ■  0  ,  R,  :  0  <  r  <  r 
1  1  1  ’  1  o 

(6l 

Lu£  +  K2ru2  "  0  »  ^2  {  ro  <  r  <  a  ’ 

where  and  are  assumed  to  be  complex  wave  numbers.  The  associ¬ 

ated  eigenfunctions  u^  and  U2  ate  assumed  to  be  complex  and  there¬ 


fore  ve  choose 
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U1 

m 

Vl  + 

lwl 

» 

u2 

m 

v2  + 

iw2 

> 

<1 

m 

nx  + 

ivl 

• 

*2 

m 

n2  + 

iv2 

» 

k2 

K1 

m 

v|) 

+  2injV 

K2 

k2 

■3 

(n 

v|)  +  2in2v 

(67) 


Substituting  into  Eq.  (66)  and  equating  real  and  imaginary  parts,  we 
get  the  four  equations 


Lvl  +  <"!-  vl^rvl  "  2ll1'01rwi 


0  , 


0  <  r  <  r 

o 

Lwx  +  (n2-  v2)rwx  +  ZrijVjTVj.  =  0  ; 

Lv2  +  (n2-  v2)rv2  -  2n2v2rw2  *  0  , 

r  <  r  <  a 
o 

Lw2  +  (n2~  \>2)rw2  +  2q2v2rv2  =  0  . 

From  these  equations  we  easily  get  the  relation 

VjLw^  -  w^Lv^  +  2n^v1r(v^  +  w£)  *  0  ,  0  <  r  <  rQ 

v2Lw2  -  w2Lv2  +  2ri2v2r(v2  +  w|)  **  0  .  vQ  <  r  <  a 

We  first  remark  that  the  operator  L  is  singular  at  r  ■  0  . 

we  consider  the  improper  Integral 


(68) 


(69) 


Therefore 
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Let  and  <2  be  the  wave  numbers  defined  by  Eqs.  (9)  and 

(11).  Aaatxmlng  £  ■  +  i^  in  these  equations,  we  get 

*2  ■  w2/c2  -  $2  „  (u2/c2  -  ^2  +  g2)  _  215^  , 

(72) 

ie|  ■  u>2 /c^  -  £2  ■  (w2/c|  -  52  +  C2)  -  2i£^2  * 

and  therefore  comparing  with  Eq.  (67)  we  find  that 

hjVi  «  n2v2  "  "  •  (73) 
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Equation  (71)  therefore  takes  the  form 


(v^  +  w^)rdr  + 


a 


(v|  + 


0  . 


(74) 


o  r 

o 

We  first  assume  that  the  quantity  inside  the  braces  is  non-zero.  Then 
Eq.  (74)  will  be  satisfied  if  either  £^  ■  0  or  m  0  .  This  says 
that  real  and  imaginary  wave  numbers  £  are  admissible.  Now  let  £ 
be  complex  so  that  £^  and  are  both  non-zero.  Then  for  the  equa¬ 

tion  to  hold  the  quantity  within  the  braces  must  be  zero.  Since  each  of 
the  integrand  is  positive-definite,  it  follows  that  u^  ■  (v^,w^)  =  0, 
and  U£  "  (v2  * W2^  2  0  *  This  implies  that  for  complex  values  of  £  , 

the  functions  u^(r)  and  u^(r)  are  identically  zero,  and  therefore 
complex  £  is  inadmissible.  We  therefore  conclude  that  frequency  equa¬ 
tion  (13)  does  not  admit  complex  wave  number  £  ,  but  can  be  real  or 

imaginary. 
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Coincidence  of  Frequencies 

It  is  easy  to  demonstrate  that  for  each  frequency  id  >  0  there 
exists  only  one  linearly  Independent  eigenfunction,  and  therefore  all 
roots  of  the  frequency  equation  are  simple  and  coincidence  cannot  take 
place.  Suppose  that  for  a  given  frequency,  Eq.  (66)  has  two  linearly 
Independent  eigenfunctions  and  $2  *  *f  we  compute  the  conjunct 

of  the  operator  and  use  the  boundary  conditions  to  evaluate  it  we  find 
that  the  conjunct  vanishes.  However,  in  this  case  the  conjunct  is  the 
Wronskian  and  vanishing  of  the  Wronskian  implies  that  the  two  solutions 
cannot  be  linearly  Independent  as  assumed.  Hence  all  the  roots  of  the 
frequency  equation  for  u  >  0  are  simple  and  coincidence  of  frequencies 


J 
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cannot  occur,  because  our  problem  is  of  Sturro-Liouville  type  [13], 

On  the  other  hand  tc  ■  0  when  u  »  £  «  0  ,  and  the  equation 
deqpnerates  to  an  Euler's  equation.  The  conjunct  of  the  operator  is 
^  ^  (d/diOWl^,^]^  "  <^»^2>  '  an*1  hence  non-zero.  We  thus  have 
a  problem  of  co-existence  when  u>  ■  5  •  0  .  In  the  theory  of  spectral 
representation  this  knowledge  is  of  fundamental  importance. 


Velocity  of  Energy  Flow 

In  an  open  connected  domain  R  ,  the  stress  equation  of  motion, 
corresponding  to  displacement  equation  (1) ,  can  be  written  as 


III e  2  3t0z  32u 

3r  +  r  Tre  +  17“  "  p 


0  , 


RtR^*  ?  , 


(75) 


where  in  terms  of  the  physical  components 


re 


r  3u  u\ 
l^3r  ~  7' 


! 

i  3u 

!  Tze  *  w  a7 


2v  ere  » 


“  2V  ez6  • 


(76) 


The- total  energy  is 


E  ■  |  e  r  dr  da  , 
R 


(77) 


where  the  scalar  energy  density  is 


*  “  ^Tre  ®re  +  T*e  Se*  +  1  pu2* 


(78) 


which  is  the  sum  of  the  potential  and  kinetic  energy  densities. 
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The  time  rate  of  change  of  energy  density  can  be  easily  obtained 
by  differentiating  the  energy  density  e  with  respect  to  the  time 
parameter  t  .  When  the  shear  stresses  satisfy  the  equation  of  motion 
(75),  the  rate  of  change  of  energy  density  is  given  by 

If  -  r]h  cirTre)  +  h  ^"ze5]  “  0  •  (79) 

Relating  the  physical  components  to  tensor  components,  we  now  define 


re  =  e  , 


urx 


re 


=  u2t 


12 


(80) 


•  ■ 

urTze  =  V  * 

-*K2 

where  u2  is  the  covariant  component  of  displacement,  t  are  the 
contravarlent  components  of  stress  tensor  density  of  weight  +1  ,  and 
e  is  the  scalar  energy  density  of  weight  +1  .  The  equation  of  con¬ 
tinuity  of  energy  and  its  flow  now  takes  a  familiar  covariant  form 


3e 

at 


VvK2> 


o  , 


K  m  1,3 


(81) 


where  we  have  used  the  standard  summation  convention,  identified  the 
coordinates  r,8,z  with  the  indices  1,2,3  respectively,  and  3^ 
represents  the  partial  derivative  with  respect  to  the  K-th  coordinate, 
[14].  On  comparing  this  equation  with  the  flow  of  matter  in  fluid,  we 
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easily  Infer  that  the  mean  velocity  of  energy  flow  Is  given  by 
t 

dt  | 

- -  .  R  :  Rx  x  R2  (82) 

jdt  |  e  d  5 
o  R 

< 

where  is  the  K-th  component  of  group  velocity. 

Details  of  the  Dispersion  Spectrum 

As  pointed  out  earlier,  the  sketching  of  the  spectrum  is  considerably 
simplified  by  the  presence  of  bounding  curves  corresponding  to  the  value 
of  parameter  Q  -*■  0+  and  Q  t  »  ,  The  spectral  lines  must  pass  through 
the  intersection  points  corresponding  to  the  family  of  bounds  Q  0+, 
and  the  family  for  Q  +  «  ,  separately,  because  the  solution  common  to  the 
two  problems  is  also  a  solution  to  the  original  problem.  In  addition,  the 
spectral  lines  cannot  cross  these  bounds  because  such  solutions  do  not 
satisfy  Eq.  (13) . 

For  convenience  of  drawing  the  spectrum  we  introduce  non-dimensional 
frequency  fi  and  non-dimensional  wave  number  X  ,  where 

•  u> /£>  ,  X  •  (a-roH/ir  ,  (83) 

and  u  ■  irc2/(a-ro)  ,  is  the  lowest  thickness-shear  frequency  of  an  infinite 
plate  of  thickness  (a-rQ)  ,  [9],  In  terms  of  non-dimensional  variables 
the  bounding  curves  for  Q  -*■  0+  are 
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n2  -  x2  -  I(l-n)/it]2  <x2  ,  q  -  0,1,2,, 


n2/c2  -  x2  -  [(i-n)/ini32  y2  »  p  -  0,1,2,, 


(84) 


where  c2  =  c2/c|  =  a/p  .  The  coordinates  of  the  Intersection  points  are 


X2 

pq 


pq 


(CRYP  "  ^2aq)^l-C2)  , 


&2<YP  '  n2“q)cR/(1_CR) 


(85) 


For  the  second  set  of  bounding  curves  corresponding  to  Q  +  ®  ,  and  the 

coordinates  of  the  Intersection  points,  we  replace  y  ,  a  by  y  ,  a  in 

p  q  m  n 

the  above  formulas.  The  values  of  y  ,  a  and  y  ,  n  for  Q  ■+■  0+  and 

p  q  m  n 

Q  t  •  ,  are  the  simple  zeros  of  Eqs.  (30)  and  (31),  respectively. 

In  the  (fl,X)-plane,  the  straight  lines  fi  -  X  for  a  -  0  ,  and 
0  ■  c^X  for  y  ■  0  ,  are  the  asymptotes  to  the  bounding  curves.  These 
asymptotic  lines  divide  the  real  plane  in  three  regions,  i)  0  >.  X  , 
ii)  X  <  fi  £.  XcR  ,  and  iii)  XcR  <  ft  <.  0  ,  where  we  have  assumed  that 
cR  <  1  .  In  the  first  region  a  and  y  are  both  real;  in  the  second 
region  a  is  pure  imaginary  and  y  is  real;  and  finally  in  the  third 
region,  a  and  y  are  both  pure  Imaginary.  However,  as  shown  earlier,  no 
roots  exist  in  the  third  region.  The  asymptotic  line  fi  «•  X  intersect 
the  bounds  Q  -*■  0+  and  Q  t  »  ,  and  their  coordinates  of  intersection 
are  given  by  Eq.  (85)  if  we  set  aQ  m  0  in  these  equations.  The  cut-off 
frequencies  are  the  roots  of  transcendental  equation  (13),  if  in 
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this  equation  we  set  o  ■  irfl/(l-n)  *  Y  ”  ryiffl/[(l-n)cR]  .  In  the  high 

frequency  region  this  equation  takes  the  simpler  form  of  Eq.  (SI) . 

■  the  dispersion  spectrum  can  now  be  calculated  using  the  above 

X 

information  as  guide.  For  real  a  and  y  ,  Eq.  (13)  must  hold.  Choosing 
appropriate  starting  values  and  increments  for  a  ,  the  corresponding  y's 
may’ be  found  by  numerical  iteration.  Each  y  can  be  Identified  to  the  mode 
it  represents  by  means  of  the  bounds  on  a  and  y  .  Once  a  and  y  are 
known,  the  corresponding  values  of  0  and  X  are  given  by 

fl2  _  (i=n)2[(yn)2  .  f 

(86) 

*2  -  (t1)^^2  -  a5]/(1/ci  - • 

As  soon  as  the  asymptotic  line  0  *  X  is  crossed,  a  becomes  imag¬ 
inary,  y  remains  real,  and  now  Eq.  (16)  must  be  satisfied.  Once  a  and 
y  are  determined  from  this  equation,  the  coordinates  0  and  X  can  be 

determined  from  Eq.  (86),  if  in  this  equation  we  replace  a  by  ia  .  As 

1 

mentioned  earlier,  when  cR  <  1  ,  there  are  no  roots  of  Eq.  (18),  when  both 
a  tend  y  are  imaginary.  Hence  the  spectral  lines  never  intersect  the 
asymptotic  line  0  •  XcR  ,  but  approach  it  asymptotically  from  above. 

Aftdr  the  solution  for  real  X  has  been  carried  out,  the  procedure  is 
repeated  for  imaginary  values  of  X  .  The  bounding  curves  are  now  circles 
and  ellipses,  rather  than  hyperbolas  as  was  the  case  for  real  X  .  The 
computational  procedure  is  the  same  as  for  real  X  ,  except  that  we 
replace  X  by  IX  in  all  relevant  formulas. 


jt 


Figures  (1,2)  illustrate  a  typical  set  of  dispersion  curves  for 
n  *  V3,  a  ■  1  and  c^  ■  ^ ,  first  in  the  (o,y) -space  and  then  in  the 
more  familiar  (Q,A)-space.  The  subscripts  1,2  on  a  and  y  indicate 
the  bounding  curves  corresponding  to  Q  +  *>  and  Q  -*■  0+,  respectively. 
The  integer  inside  the  brackets  indicates  the  'mode'  or  value  of 
p,q,m  or  n,  for  example  a2(3)  =  aq=3*  Yl^  =  Yn=l 


,  etc. 
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Figure  Captions 


1  Dispersion  Spectrum  in  ( a  , y  )  space  with  Bounds  and 
Cutoff  Frequencies. 

2  Dispersion  Spectrum  in  (  A  ,  fi  )  space  for  Real  and  A 
Imaginary  A  with  Bounds  and  Asymptotes. 


Id 


if] 


SKCURITY  CLASSIFICATION  OF  THIS  PAGE  fWTwri. 


REPORT  DOCUMcNTATiCM  PAGE 

KKAO  JNSTkUCTlONS 
ilKl  Ok’I-:  COMF1, l-TlNf.  KOi?M 

1.  REPORT  NUMBER 

107 

mmrn 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (atuf  Subtitle) 

5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

TORSIONAL  WAVES  IN  A  BIMETALLIC  CYLINDER— 

PART  1  •  AYTAT.T.Y  HflMnORNFflTTS  P.ASE 

TECHNICAL 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^ 

0.  CONTRACT  OR  GRANT  NUMBER^; 

R.K.  KAUL,  R.P .  SHAW,  W.  MULLER 

N00014-71-C-0108  * - 

a.  PERFORMING  organization  name  and  aooress 

DEPT.  OF  ENG.  SCI.,  AERO.  ENG.  AND  NUC.  ENG. 

SUNY  AT  BUFFALO,  BUFFALO,  NY  14214 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

12.  REPORT  DATE 

OCTOBER  1979  - - 

OFFICE  OF  NAVAL  RESEARCH  -  STRUCTURAL  MECHANICS 

13.  NUMBER  OF  PAGES 

38  PAGES  +  2  FIGURES 

14.  MONITORING  AGENCY  NAME  6  ADDRESS^!/  dlllcront  from  Controlling  Olllcc ) 

15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

15*.  OECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 

to.  DISTRIBUTION  STATEMENT  (of  this  Report) 

APPROVED  FOR  PUBLIC  RELEASE.  DISTRIBUTION  IS  UNLIMITED. 

17.  DISTRIBUTION  STATEMENT  f of  tba  abstract  entorad  fn  Block  30.  II  dllloront  from  RaportJ 

la.  SUPPLEMENTARY  NOT£S 

19.  KEY  WORDS  (Continue  on  reverse  aide  it  necessary  and  Identity  by  block  number) 

TORSIONAL  WAVES;  COMPOSITE  STRUCTURE 

20.  ABSTRACT  (Continue  on  reverse  aide  If  necesesry  and  Identify  by  block  number) 

The  dispersion  spectrum  is  found  for  axially  symmetric  torsional  waves  in  an 
elastic,  bimetallic  rod  with  cylindrical  core  and  concentric  outer  casing. 
Plotting  of  the  various  branches  of  the  spectrum  is  simplified  by  the  pre¬ 
sence  of  discrete  invariant  points  which  are  independent  of  the  material 
properties,  and  through  which  the  spectral  lines  must  pass.  The  slope  and 
curvature  of  the  spectral  lines  at  cut-off  frequencies,  the  asymptotic 
approximations  at  high  frequency,  the  non-existence  of  complex  branches,  the 
problem  of  co-existence  and  the.. concent,  of  energy,  flow. are  .also  studied..,. 

DD  i  j an^73  14/3  EOITION  O F  1  NOV  6&  IS  OUSOl.ETE 


S  N  0102-  l.r-  0  1  .t-  6601 


SECURITY  CL  ASSIHCAYlOItToF  Th7s  PMiK  Oat*  Snlctmd) 


